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Layered organic superconductors of the BEDT family are model systems for understanding the interplay of 
the Mott transition with superconductivity, magnetic order and frustration, ingredients that are essential to un¬ 
derstand superconductivity also in the cuprate high-temperature superconductors. Recent experimental studies 
on a hole-doped version of the organic compounds reveals an enhancement of superconductivity and a rapid 
crossover between two different conducting phases above the superconducting dome. One of these phases is a 
Fermi liquid, the other not. Using plaquette cellular dynamical mean field theory with state of the art continuous¬ 
time quantum Monte Carlo calculations, we study this problem with the two-dimensional Hubbard model on the 
anisotropic triangular lattice. Phase diagrams as a function of temperature T and interaction strength U/t are 
obtained for anisotropy parameters t' = 0.4t, t' = O.St and for various fillings. As in the case of the cuprates, 
we find, at finite doping, a first-order transition between two normal-state phases. One of theses phases has a 
pseudogap while the other does not. At temperatures above the critical point of the first-order transition, there 
is a Widom line where crossovers occur. The maximum (optimal) superconducting critical temperature TJ^ at 
finite doping is enhanced by about 25% compared with its maximum at half-filling and the range of U/t where 
superconductivity appears is greatly extended. These results are in broad agreement with experiment. Also, 
increasing frustration (larger t'/t) significantly reduces magnetic ordering, as expected. This suggests that for 
compounds with intermediate to high frustration, very light-doping should reveal the infiuence of the first-order 
transition and associated crossovers. These crossovers could possibly be even visible in the superconducting 
phase through subtle signatures. We also predict that destroying the superconducting phase by a magnetic field 
should reveal the first-order transition between metal and pseudogap. Finally, we predict that electron-doping 
should also lead to an increased range of U/t for superconductivity but with a reduced maximum Tc. This 
work also clearly shows that the superconducting dome in organic superconductors is tied to the Mott transition 
and its continuation as a transition separating pseudogap phase from correlated metal in doped compounds, 
as in the cuprates. Contrary to heavy fermions for example, the maximum Tc is definitely not attached to an 
antiferromagnetic quantum critical point. That can also be verified experimentally. 

PACS numbers: 74.70.Kn, 71.10.Fd, 71.30.+h, 74.20.Mn 


I. INTRODUCTION 

In organic charge transfer salts, such as 
/^-(BEDT-TTF) 2 X (/^-(ET) 2 X ) or EtnMe4-nPn[Pd(dmit)2]2 
(Pd(dmit) 2 ), a first-order phase transition between a super¬ 
conductor and a Mott insulator is induced by pressure [1-5]. 
These materials also present a wide range of intriguing phe¬ 
nomena such as unconventional superconductivity, magnetic 
ordering [5], pseudogap [6] valence-bond solid phases [3] 
and some of them are even spin-liquid candidates. [7, 8]. 
Excellent reviews are available. [9, 10] 

The presence of the Mott transition and of spin-liquid states 
in the phase diagram suggests that strong electronic correla¬ 
tions and electronic frustration are key to the physics of the 
organics. The one-band Hubbard model on an anisotropic tri¬ 
angular lattice near half-filling is the simplest model that cap¬ 
tures this physics [9-1 1], although consensus has not yet been 
completely reached [12, 13] on this point. 

Unraveling the physics of these layered materials should 
also be helpful to shed light on cuprate high-temperature 
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superconductors. Indeed, these two classes of materials 
give a complementary perspective on the crucial role of 
the Mott transition. In the organics, the Mott transition is 
bandwidth-controlled whereas it is doping-controlled in the 
cuprates. [14, 15] The analogy between these two classes 
of materials has been reinforced through recent experimental 
studies [16] of doped organics [17, 18] that show a supercon¬ 
ducting dome [19] as a function of pressure as well as a rapid 
change from non-Eermi liquid (pseudogap) to Eermi-liquid 
like metal at a critical pressure in the normal state. [16, 20] 
These analogies motivate our work. In short, our calculations 
explain these different features just as calculations performed 
with the same methods explain many of the key features of the 
cuprate phase diagram. [21] In particular a first-order transi¬ 
tion at finite doping and its associated Widom line play a cru¬ 
cial role as in the cuprates. [22, 23] We define superconduc¬ 
tivity as strongly-correlated when it arises in the presence of 
interactions larger than, or of the order of, those necessary to 
lead to a Mott transition at half-filling. 

Some of the striking experimental results that we address 
and understand theoretically in this paper come from re¬ 
cent work on /^-(ET) 4 Hg 2 . 89 Br 8 by Oike et al [16]. That 
compound is considered as a 10% doped analog of n- 
(ET) 2 Cu 2 (CN )3 . The main observation is that maximally 
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enhanced superconductivity and a normal-state crossover to a 
non-Fermi liquid phase appear concomitantly around a pres¬ 
sure where mobile carriers decrease rapidly. Also, the range 
of pressures spanned by the superconducting dome in k- 
(ET) 4 Hg 2 . 89 Br 8 is about six times the range where it appears 
in the half-filled analog /^-(ET) 2 Cu 2 (CN )3 . 

We work with the Hubbard model on the anisotropic 
triangular lattice as a function of temperature, interaction 
strength U/t and filling n, for different values of frustration 
characterized by the ratio of near-neighbor hopping t'/tin dif¬ 
ferent directions. Values for t'/t are inspired from Kandpal et 
al. [24] and from Nakamura et al. [25] who found, using ah 
initio density functional theory, that /^-(ET) 2 Cu[N(CN) 2 ]Cl 
could be modeled by t' = 0.4t and /^-(ET) 2 Cu 2 (CN )3 by 
t' = O.St. Comparing the results for two values of t'/t helps 
understand the effect of magnetic frustration on the phase di¬ 
agram. 

We use a cluster generalization of dynamical mean-field 
theory. [26-29] This approach has already led to numerous 
results that can be confronted with experiments, thus permit¬ 
ting to address important issues in cuprates, Pd(dmit) 2 and n- 
(ET) 2 X such as the pseudogap, superconductivity, Mott tran¬ 
sition, magnetic ordering, thermodynamic properties, and un¬ 
usual criticality in organic compounds. [15] The assumption 
inherent in our approach is that the main physics of the or¬ 
ganics originates from strong correlations that occur at short 
distances due to on-site repulsion U and near-neighbor su¬ 
perexchange J that are present in the Hubbard model. This 
is the assumption behind the state of the art method that we 
use for a 2 X 2 cluster embedded in a self-consistent dynam¬ 
ical mean-field. Larger cluster calculations would be neces¬ 
sary if this assumption was proven incorrect. The agreement 
between our results at low temperature and those obtained re¬ 
cently through variational Monte Carlo methods helps estab¬ 
lish the validity of the approach. [30] Agreement with several 
experimental facts strongly suggests that our approach is rele¬ 
vant for experiment. We make predictions for experiment that 
can falsify the theory. In the absence of an exact solution to 
the Hubbard model, we are making a minimum number of as¬ 
sumptions and suggesting experiments that can falsify them 
by disagreeing with the consequences. 

The link between the normal state of high-temperature 
cuprate superconductors and that of the organics is illustrated 
schematically [31] in Eig. I. Disregarding temperature, the rel¬ 
evant variables are interaction strength U/t , doping J, and 
frustration as measured by the ratio of second to first nearest- 
neighbor hopping t'/t. Cuprate superconductors, in the red re¬ 
gion, are easily doped but are little influenced by pressure and 
the range of t'/t varies little between different compounds. 
Layered organic compounds on the other hand are half-filled, 
with a broad range of possible values of frustration t'/t from 
compound to compound, and their bandwidth to interaction 
ratio is strongly influenced by pressure. A pseudogap appears 
through a second-order transition when doping is increased 
from the yellow region containing the Mott insulating phase 
at half-filling. The pseudogap phase ends in a first-order man¬ 
ner on the magenta surface in Eig. I. The latter first order tran¬ 
sition extends from the first-order Mott insulator-metal tran- 



FIG. 1. Schematic generalized normal state phase diagram for the 
cuprates and the layered organics, in the limit of low temperature 
neglecting broken-symmetry phases. The yellow surface represents 
the Mott insulating phase at zero hole-doping (5 = 0, or half-filling 
n = 1. Leaving the Mott insulating phase by increasing doping 
yields a second-order transition to a pseudogap phase up to the ma¬ 
genta surface where a first-order transition to a correlated metal oc¬ 
curs (at least close to half-filling). That first-order transition extends 
from the first-order Mott metal-insulator transition that occurs at the 
boundary of the yellow region. Cuprate superconductors are found 
in the red region at negative t'/t. The half-filled layered organics 
are found for different t' /t along the blue regions in the zero hole¬ 
doping plane. They are strongly influenced by pressure, whose ef¬ 
fect is represented vertically, although the ratio t' /t will generally 
also be influenced by pressure. A doped layered organic is found in 
the green region. It allows one to expand the analogy between or¬ 
ganics and cuprates. The organic and cuprate lattices are different, 
as illustrated respectively by the left and right insets. Nevertheless, 
the physics of interactions, frustration and doping is present in both 
types of compounds. Superconductivity and antiferromagnetism are 
broken symmetry phases that are strongly influenced by the underly¬ 
ing normal state illustrated by this flgure. In particular, the supercon¬ 
ducting phase has a maximum Tc in the vicinity of the boundary of 
the magenta region. 

sition occurring at half-filling at the boundary of the yellow 
region. 

The range of parameters that could be relevant for the 
doped organic compound /^-(ET) 4 Hg 2 . 89 Br 8 is indicated by 
the green region in Fig. I. This compound then offers the in¬ 
teresting possibility to investigate the pseudogap to correlated 
metal transition in the normal state by cutting the first-order 
magenta surface along a direction different from that of the 
cuprates. The effect of this transition on the superconducting 
state is also a key question that we address here. 

The Hubbard model and the Cellular Dynamical Mean- 
Field Theory on a plaquette with Continuous-Time Quantum- 
Monte-Carlo impurity solver are presented in Sec. II. This 
work would not have been possible without recent improve¬ 
ments of this solver related to sign problem minimization, [32] 
ergodicity [33], and speedup. [34] We begin in Sec. Ill with 
a short summary of some previous results for the cuprates. 
We then present in Sec. IV results for the normal state, show¬ 
ing the Widom line that emerges from the first-order transi¬ 
tion [35, 36] between a pseudogap phase and a metal. [22] 
This plays a crucial role for the cuprates. [21] The results 
in Sec. V are for two different lattice anisotropies, t'/t, or 
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equivalently, frustration. We investigate the Neel antiferro¬ 
magnetism (AFM) and d-wave superconductivity (SC) on the 
same footing but the relative stability of the phases is not stud¬ 
ied. For t' jt = 0.4, half-filling, 1%, 10% hole-doping and 
10% electron-doping are investigated. We find for a 1% hole¬ 
doping that the maximum (optimal) superconducting critical 
temperature (TJ^ ) as a function of pressure or (interaction 
strength) is enhanced by approximately 25% and the range of 
superconductivity is multiplied by a factor of six on the pres¬ 
sure axis it/U). The range of pressure where superconductiv¬ 
ity exists for 10% doping is similar. We also obtain the 
line in the T — f/ — n phase diagram. The case t'/t = 0.8 
is considered only for the 10% hole-doped case and at half¬ 
filling due to a worse sign problem. Discussions in Sec. VI 
include the role of long wave-length fiuctuations on broken- 
symmetry phases, the role of antiferromagnetic quantum crit¬ 
ical points and of the Mott transition on the superconducting 
dome, contact with experiment, predictions, limitations of the 
approach, and perspectives. The most important conclusions 
are summarized in Sec. VII. 



FIG. 2. Periodic partitioning of the anisotropic triangular lattice into 
2x2 frustrated square clusters for this work using CDMFT. 


II. MODEL AND METHOD 

We consider the single band Hubbard Hamiltonian on the 
anisotropic triangular lattice in two dimensions 

H tij 4^Cj^ + ^ X] X] 

i i,(J 

where Uj is the hopping amplitude between neighboring sites, 
Cia- and cj^ respectively destroy and create an electron of spin 
a at site i, riicj is the density of electrons of spin a at site i, 
/i is the chemical potential and U is the on-site Coulomb re¬ 
pulsion. This model was proposed several decades ago for the 
organics. [11] Its validity has been revisited recently for spe¬ 
cific compounds. [37] In particular, it has been argued that 
near-neighbor repulsion V was important. Since previous 
studies for the cuprates have shown that V does not infiuence 
the phase diagram in a dramatic way in the strong-correlation 
limit, [38] we neglect this term in this initial study. 


As illustrated in Fig. 2, we take tij = t for nearest-neighbor 
bonds, and Uj = t' for the diagonal bond. The isotropic trian¬ 
gle is recovered for t' = t. The ratio t'/t is both an anisotropy 
parameter and a measure of magnetic frustration. The two 
expressions are used interchangeably. In the figures, inverse 
temperature P is given in units of 1/t. 

Tc stands for the superconducting critical temperature. It 
was called in Ref. 39 to emphasize that it is the dynamical 
mean-field transition temperature, that differs from the true 
superconducting transition temperature. We use for the 
maximum value that this quantity takes at a given doping as a 
function of pressure or interaction strength. 

A. CDMFT 

In two dimensions, momentum-dependence of the self¬ 
energy is important. Cellular dynamical mean field the¬ 
ory [40] (CDMFT) for the Hubbard model takes into account 
short-range correlations in addition to interaction-induced dy¬ 
namical correlations; single-site DMFT [26] is not appropriate 
to study the momentum-dependence associated to d-wave su¬ 
perconductivity. The key approximation is to restrict the self¬ 
energy to a local cluster and neglect its spatial dependence 
beyond the cluster [28, 29, 41, 42]. 

In practice, CDMFT embeds a cluster of finite size in a 
non-interacting electronic bath. The impurity problem (clus¬ 
ter and bath) is then solved and the bath is determined self- 
consistently by demanding that the lattice Green function pro¬ 
jected on the cluster equal the Green function obtained from 
the impurity problem. 

To be more specific, the lattice Green function in Matsubara 
frequencies is obtained from 

k) = 3“ /^) I ^(k) ^clip^n) (2') 

where k is the wave vector associated with translational in¬ 
variance from cluster to cluster, fi is the chemical potential, 

f(k) the full hopping matrix (including intra-cluster and inter¬ 
cluster hoppings) and is the self-energy of the clus¬ 

ter, imposed to be equal to the self-energy of the lattice in the 
CDMFT approximation. The hat on symbols specifies that 
they are matrices in the basis of cluster Wannier states. 

The self-energy is related to the Green function of the clus¬ 
ter through Dyson’s equation 

- Sd (3) 

where the free propagator on the cluster is defined by 

Qq (iujji) = {iuJn H“ /i) i ^ (4) 

with the one-body part of the Hamiltonian and Acz the 
hybridization function that defines the bath and its coupling to 
the cluster. 

The projection of the lattice Green function on the clus¬ 
ter leads to the self-consistent equation for the hybridization 
function 

/ dk - ~ 

(27r)2 k), (5) 
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where is the number of sites on the cluster and the in¬ 
tegral is over the reduced Brillouin zone. A high frequency 
expansion of both sides of this equation proves that is the 
hopping matrix within the cluster. [43] 

The uniform spin susceptibility, or Knight shift, is defined 
by 

Sz = l (A^t - N^) 

r0 

Xz{q = 0,UJ = 0) = / {Sz{t)Sz( 0)) dr, (6) 

Jo 

with and the total number of up and down spins, re¬ 
spectively, on the cluster. 

The calculation starts from a guess for the hybridization 
function Ad . This gives the so-called dynamical mean-field, 
Eq. 4, namely the free propagator for the cluster. Assum¬ 
ing that the impurity problem can be solved, Dyson’s Eq. 3 
then gives the self-energy which is needed to obtain the lat¬ 
tice Green function, Eq. 2, entering the right-hand side of the 
self-consistency Eq. 5. Using Dyson’s equation again on the 
left-hand side of that self-consistency equation leads to a new 
guess for the hybridization function. This process is iterated 
until convergence. 

Calculating the cluster impurity Green function Gd is the 
difficult problem. This is done here using a continuous-time 
quantum Monte Carlo method (CTQMC). 


are both block diagonal. The largest block is 2 x 2 because the 
Ai representation occurs twice. The calculations presented 
here are possible only if the angle defining rotations in this 
2x2 block is chosen to minimize the sign problem. [32] In 
addition it is necessary to use a modification of the original 
algorithms to ensure ergodicity in the presence of d-wave su¬ 
perconductivity. [33] We also speedup the calculation with the 
Lazy Skip List algorithm. [34] 

In normal phase studies, this basis respects the symmetries 
of the lattice that are compatible with the partitioning. In bro¬ 
ken symmetry phases, such as magnetically ordered or d-wave 
SC, symmetry-breaking is allowed only for the hybridization 
function. The cluster continues to respect the original Hamil¬ 
tonian symmetries. There is no mean-field factorization on the 
cluster. 

Other popular continuous-time quantum Monte Carlo im¬ 
purity solvers involve expansion in powers of the interac¬ 
tion. They have better scaling than CT-HYB with cluster size. 
However, they need very large order expansion at large U/t, 
which makes them converge slowly, and they have a severe 
sign problem at large interaction strengths U/t and frustration 
f/t. [45] 


C. Broken symmetry phases 

The Green function for superconductivity is written in 
Nambu notation as 


B. Hybridization-Expansion Continuous-Time Quantum 
Monte Carlo 


CTQMC provides a statistically exact solution of the impu¬ 
rity problem exempt from imaginary-time discretization error. 
The large values of U/t , low temperatures and large frus¬ 
tration that we need can be attained only with the hybridiza¬ 
tion expansion algorithm (CT-HYB). [44] Extensive reviews 
of CTQMC solvers are available. [45-47] 

The cluster that tiles the infinite anisotropic triangular lat¬ 
tice is illustrated in Eig. 2. It allows a singlet ground state. 
To speedup the calculations, one chooses a single-particle ba¬ 
sis that transforms as the irreducible representations of the 
cluster-Hamiltonian symmetries. [47] The point group sym¬ 
metry C 2 v of tho anisotropic cluster as well as charge and spin 
conservation lead to the following single-particle basis 


^Aia — (^Icr 3“ 

C'Aia ~ (^2cr + Qcr) 

(^Icr ^3cr) 

^B2cr — (^2cr C4cr) , 


(7) 


where the indices are those of Eig. 2 and where Ai, Bi, B 2 
are irreducible representations of C 2 v^ ^2 being empty. In this 
basis, the hybridization function A and cluster Green function 


. ( 8 ) 

with = (c|, C|) where c| and are row vectors as de¬ 
fined by equation Eq. 7. The d-wave superconducting order 
parameter transforms as the A 2 representation of the C 2 V sym¬ 
metry group. Hence only entries in the Gork’ov function ^(r) 
transforming as A 2 can be finite, e.g. for singlet pairing 

(^) : = - < TrCBif(T)cB2l > 

= - < TtCB2^{-'^)cBiI > (9) 

=: J^B2,B^{-r) 

To determine the region where the SC phase is allowed we 
calculate the order parameter 

dsc •= (10) 

Since antiferromagnetism does not break C 2 V symmetry, no 
additional entry is needed in the Green function matrix. We 
only need to let up and down spins take independent values. 


III. A BRIEF REVIEW OF QUANTUM CLUSTER 
RESULTS FOR THE CUPRATES 

In this section we only briefiy summarize some of the main 
results obtained with cluster generalizations of DMET for the 
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cuprates and give a few representative references. A more 
detailed but not exhaustive review can be found in Ref. 15. 

There are two cluster generalizations [28, 29, 41] of DMFT. 
We described CDMFT above. In the Dynamical Cluster Ap¬ 
proximation (DCA) [48] the clusters are built in momentum 
space. Whatever the method used, all groups have found 
a pseudogap in the normal state near half-filling. [22, 49- 
54] The Mott transition at half-filling is first-order, [55, 56] 
and large cluster studies find [27] a d-wave superconduct¬ 
ing transition temperature Tc at finite doping. The zero- 
temperature order parameter has a dome shape instead of in¬ 
creasing monotonously towards half-filling when the interac¬ 
tion strength is large enough that there is a Mott insulator at 
half-filling or when antiferromagnetism is allowed to compete 
with superconductivity. [57-63] The larger cluster studies [62] 
find that the dome ends at a finite doping away from half¬ 
filling. 

All methods agree with the existence of crossovers at high- 
temperature associated with the opening of a pseudogap as 
doping is reduced towards half-filling. On 2 x 2 clusters with 
the CT-HYB solver it was possible to reach lower tempera¬ 
tures than previous studies. Scanning chemical potential over 
a very fine mesh allowed the discovery of a first-order transi¬ 
tion in the normal state at finite doping. [22] That first-order 
transition ends at a critical point that is continued as a Widom 
line in the supercritical region. That first-order transition with 
its Widom line then becomes an organizing principle for the 
observed crossovers and for the superconducting dome. The 
Widom line is described in more details in the following sec¬ 
tion that begins our discussion of organics. Earlier work on 
normal state first-order transitions at finite-doping is discussed 
at the end of that section. 


IV. NORMAL STATE PSEUDOGAP, FIRST-ORDER 
TRANSITION, AND WIDOM LINE 

The first-order Mott transition at half-filling is well docu¬ 
mented. [64, 65] The blue shaded region in Fig. 5(a) identifies 
the region of parameter space where normal-state hysteresis is 
found. Metallic and insulating states there can coexist. The re¬ 
sults are similar to those obtained at half-filling in the unfrus¬ 
trated case t' = 0. [55, 56] The positive slope of the transition 
in the T — (7 plane (negative in the T — t/f/ plane) comes 
from the smaller entropy of the insulating phase compared 
with the metallic phase, as deduced from Clausius-Clapeyron 
arguments. [1, 36] Indeed, from dG = —SdT + DdU + iidn, 
where G is the Gibbs free energy, S the entropy, D double 
occupancy, and n the filling, we find that the slope of the tran¬ 
sition line is 


dT Di — Dm 
W ^ Si-Sm 


( 11 ) 


where index / is for the insulating phase and M for the metal¬ 
lic one. The smaller entropy of the insulating phase comes 
from the tendency to form local singlets. [1, 22, 35, 36] 

In this section, we focus on the less familiar first-order tran¬ 
sition found at finite doping. Consider the case t' = OAt 


U/t 



FIG. 3. Double occupancy Docc as a function of pressure (bottom 
horizontal axis) or interaction strength U/t (top horizontal axis) for 
fixed filling, n = 0.99. The value t = 0.044 eV is used to convert to 
physical units. [64] The lower horizontal axis is labeled t/U io sug¬ 
gest the pressure dependence, but the numbers on that horizontal axis 
are given by the value of 1/U expressed in electron-Volt using the 
above conversion factor. At T = t/60 there is a first-order hysteresis 
region: the brown squares are obtained from the insulating solution 
and the blue squares for the conducting solution. At T = t/12, there 
is no hysteresis, only an inflexion point that determines the Widom 
line. 



FIG. 4. The imaginary part of the Matsubara Green func¬ 
tion Im. G{iun) plotted as a function of Matsubara frequency 
gives information about the density of states at the Fermi level, 
—21m. G{iuJn ir]) with 77 ^ 0. The behavior differs depend¬ 
ing on the value of U/t . The decrease towards zero or the density 
of states at larger U/t indicates a pseudogap, while its increase at 
smaller U/t indicates a metallic phase. 


at 1% doping, namely n = 0.99. Fig. 3 shows two jumps 
of double-occupancy delimiting a coexistence region at low 
temperature, and a smooth dependence on U/t at high tem¬ 
perature. The jumps at low-temperature define the hysteresis 
region of a first-order transition. The low pressure phase ex¬ 
hibits a pseudogap while the high-pressure phase is a more 
standard metal. The inflexion point at high temperature de¬ 
fines a crossover. The locus of these infiexion points is associ¬ 
ated with the so-called Widom line of the first-order transition. 
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In the theory of fluids, the Widom line is defined as the line 
where the maxima of different thermodynamic response func¬ 
tions touch each other asymptotically as one approaches the 
critical point of the first-order transition. [66] Investigations 
on the phase diagram of fluids have shown drastic changes 
in the dynamics upon crossing the Widom line. [66, 67] By 
analogy, in the cuprates the Widom line has been identified as 
the organizing principle for the pseudogap and resulting phase 
diagram of the high-temperature superconductors. [21-23, 68] 

In the present context, along the Widom line a crossover 
from a metallic state to a pseudogap metal is also seen. This is 
illustrated by the frequency dependence of the imaginary part 
of the local Matsubara Green function in Fig. 4 

at T slightly above the critical point of the first-order tran¬ 
sition. For large values ofU/t , lm{G{iuJri)) aims upwards 
as frequency decreases, indicating a small density of states 
at the Fermi level, consistent with a pseudogap. By contrast, 
for smaller values of U/t, Im(G(icc;n)) extrapolates to a finite 
density of states at the Fermi level, consistent with an ordinary 
metal. We suggest that this crossover corresponds to the one 
seen experimentally in doped organics. [16, 20] At lower tem¬ 
perature, the transition between the pseudogap metal and the 
more ordinary metallic phase occurs discontinuously through 
the first-order transition illustrated in Fig. 3. 

Fig. 5(b) displays the normal state phase diagram at 1% 
doping. There is a coexistence region, in blue, coming from 
the first-order transition and a Widom line that extends above 
the critical point of that first-order transition. There is also a 
Widom line in the half-filled case (not shown). In the context 
of the cuprates, this first-order transition was found at fixed 
f//t as a function of doping. [22, 35, 36] The results in those 
papers [22, 35, 36] clearly show a surface of first-order tran¬ 
sitions that is continuously connected to the Mott transition 
at half-filling. [69] The results of Fig. 5(b) are in a way a 
constant doping cut of the finite t' version of that first-order 
surface. The critical point in Fig. 5(b) occurs at a tempera¬ 
ture about 60% lower than the corresponding temperature at 
n = 1. That rapid drop is also observed in the square-lattice 
results. 

We end this section with brief comments on early work on 
the doping-induced Mott transition of the Hubbard model. In 
single-site DMFT, [70] it was found for t' = 0 that upon dop¬ 
ing there is a first order transition between a half-filled Mott 
insulator and a finite-doping metal. Essentially the same result 
was found with DC A for the square lattice and various positive 
t'/t (electron-doped case in the language of cuprates). [71] 
Later work with the same methods [72] suggested that at 
t' = 0 there is a quantum-critical point instead of a first-order 
transition but the lowest temperature reached was large com¬ 
pared to those where the first-order transition was found in 
CDMFT with CTQMC solver. [22, 35, 36] CDMFT studies 
with an exact diagonalization solver [52] also found a clear 
first-order transition for positive t'/ton the cuprate square lat¬ 
tice. As in Refs. [22, 35, 36], that first-order transition sepa¬ 
rates a pseudogap phase and a metal instead of a separating a 
Mott insulator and a metal as found in the above early DCA 
study. [71] In that study, [71] it was noted that the phase tran¬ 
sition appeared only for U/t larger than the bandwith. How¬ 


ever, it was in the work of Refs. 22, 35, and 36 that the critical 
end line of the finite-doping surface of first-order transitions in 
the {U,T,S) space of parameters was shown to be connected 
to the critical end point of the Mott transition at half-filling. 
This is an important step to differentiate strong and weak cor¬ 
relation effects. [15] The pseudogap phase appears at finite 
doping only if there is a Mott insulator at half-filling. At zero 
temperature the transition between the Mott insulator and the 
pseudogap phase is second order. [35, 36] The pseudogap is 
then different from the Mott gap even though they both ap¬ 
pear in the generalized phase diagram. [22] The significance 
and existence of the Widom line was noted in Refs. 22, 23, and 
68 . 


V. PHASE DIAGRAM, INCLUDING 
SUPERCONDUCTIVITY AND ANTIEERROMAGNETISM 

We present phase diagrams for the normal state, Neel an¬ 
tiferromagnetism (AFM) and d-wave superconductivity (SC). 
In CDMFT, AFM and SC fiuctuations are treated on equal 
footing both on the cluster and in the bath but we only allowed 
one symmetry-breaking at a time. Values for t'/t are inspired 
from Kandpal and al. [24] and from Nakamura and al. [25] 
Using ab initio density functional theory, they found that n- 
(ET)2 Cu[N(CN)2 ]C1 could be modeled by t' = 0.4t and n- 
(ET) 2 Cu 2 (CN )3 by t' = 0.8t. However, extensive Huckel 
calculations had previously found higher frustration (t'/t) val¬ 
ues for these compounds [73]. We have to keep in mind these 
uncertainties when we make contact with real materials. 

We present first the case t' = 0.4t, then t' = 0.8t. The 
intermediate frustration t' = OAt results are presented first 
because the sign problem is less severe in that case, allow¬ 
ing a more thorough investigation of the SC phase diagram. 
Even though the SC phase would be mostly hidden by antifer¬ 
romagnetism in this case, we find that the results for the pure 
SC phase (forbidding AFM) are qualitatively similar to the re¬ 
sults we discuss for larger frustration t' = O.St. In the later 
case, AFM is generally sufficiently suppressed that its neglect 
is justified. We are restricted to commensurate antiferromag¬ 
netism. We also did not allow non-collinear spin order. A 
more in-depth discussion of magnetic order and its impact on 
the presence of SC in real compounds can be found in subsec¬ 
tions V B and VI A. 


A. t' = 0.4t 

Our results for the phase diagrams at different dopings 
are summarized in Fig. 5. A few more properties for the SC 
phase are displayed in Fig. 8. Let us discuss the various phases 
in turn. 

The phase diagram at half-filling Fig. 5(a) depicts the 
phase diagram at half-filling. The blue region delimits the 
metal-insulator coexistence region associated with the first- 
order Mott transition in the normal phase. The critical value 
ofU/t agrees within error bars with previous results. [12, 74] 
The d-wave SC phase is observed in proximity to the normal- 


7 



(a) (b) 


U/t u/t 

35. 17. 11. 8.4 6.7 5.6 35. 17. 11. 8.4 6.7 5.6 



(c) (d) 

FIG. 5. (color online) Phase diagrams for the Hubbard model on the anisotropic triangular lattice with t' = 0.4t for various fillings. The 
antiferromagnetic state has been studied for the half-filled case and for 10% doping. The value t — 0.044 eV is used to convert to physical 
units. [64] The lower horizontal axis is labeled t/U io suggest the pressure dependence, but the numbers on that horizontal axis are given by 
the value ofl/U expressed in electron-Volt using the above conversion factor. The same convention is used throughout the paper. Lines are 
guides to the eye. 

(a) Phase diagram for the half-filled case. In the blue region, the Mott insulator and the metallic state coexist. At the blue points, we find a first- 
order jump in the normal state double occupancy. The lines between the points can be identified as the spinodal lines, d-wave superconductivity 
occurs in the orange region: Decreasing pressure, we find a first-order jump of the SC order parameter to zero at the orange points while upon 
increasing pressure we find a second-order transition. The AFM phase occurs below the red line that interpolates between the triangles where 
we find the Neel second-order transition. 

(b) Phase diagram for the 1% hole-doped case. The colors have the same meaning as for the half-filled case, except that in the blue region 
two different metallic states are found instead of a metallic and an insulating state like at half-filling. First-order jumps are observed at the 
blue points. There is coexistence in the blue region. Also, the red dots connected to the blue region indicate a strong crossover (Widom line) 
between a pseudogap state at small pressure (large U/t) and a metallic state at large pressure (small U/t). Orange points denote where we 
detect a second-order transition from the SC state to the normal state. 

(c) Phase diagram for the 10% hole-doped case. The transitions between the normal and SC state (orange circles) are second-order. The AFM 
phase is between the red lines. These lines interpolate between the second-order Neel transition that we find where the triangles are located. 

(d) Phase diagram for the 10% electron-doped Hubbard model. The brown dot and line are extrapolations. The transitions between the normal 
and SC state (orange circles) are second-order. 
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state first-order Mott transition. When pressure is increased in 
the SC phase, (interaction strength is decrease), it disappears 
in a second-order manner. The zero-temperature results ob¬ 
tained previously [74] suggest that in that limit the SC phase 
in Fig. 5(a) will extend beyond the phase boundary for AFM. 
When pressure is decreased, SC gives way to the insulating 
phase through a first-order jump. The maximum SC critical 
temperature (TJ^ ) is attained close to the Mott transition to 
the insulator. All the qualitative results agree with previous 
theoretical studies on the unfrustrated square lattice at finite 
temperature [39] as well as with experimental observations in 
various organic compounds of the k-(ET) 2 X family and of the 
Pd(dmit) 2 family [1-3]. 

Note that the low-pressure boundary (spinodal line) where 
the metastable metallic phase disappears discontinuously in 
favor of a stable insulator does not coincide with the low 
pressure boundary where the SC phase disappears discontinu¬ 
ously. The two boundaries are however in very close proxim¬ 
ity. There is no reason for the regions of metastability of the 
normal and superconducting phases to exactly coincide. 

When AFM order is permitted, it dominates a wide area in 
the temperature-pressure (T — P) plane. The maximum of 
the SC dome does not coincide with an AFM quantum critical 
point, as can be seen for example in 5(c). Further results on 
AFM appear in Fig. 9(a) and are discussed in Secs. VB and 
VIA. Zero-temperature studies obtained with CDMFT [74] 
and with other methods [75, 76] do suggest that for this value 
of f/t AFM is the most stable magnetic phase. 

Phase transition between pseudogap phase and metallic 
phase For small hole doping (1% in Fig. 5(b)) one finds 
the first-order transition discussed in Sec. IV. The region 
where hysteresis is found is indicated by the blue region in 
Fig. 5(b). Comparison with half-filling in Fig. 5(a) demon¬ 
strates that this transition is continuously connected to the 
first-order metal-insulator Mott transition and that it occurs at 
larger interaction strength as doping is increased, as suggested 
by the border of the magenta region in Fig. I. In other words, 
upon doping, the Mott insulator evolves continuously into a 
conducting pseudogap phase different from the metallic phase 
at larger pressure. This is in complete analogy with the results 
found previously [35, 36] for the unfrustrated square lattice. 
[77] Our results are qualitatively similar to those of Fig. 1 of 
Ref. 35 for the square lattice. We indeed also find that the 
first-order transition occurs at lower temperatures as doping 
is increased and is not accessible to our simulations for hole- 
dopings as small as 4% for temperatures down to T/f = 1/60 
and interaction strength up to P/t = 20. 

The type of pseudogap discussed here is a strong- 
correlation effect as follows from the fact that it appears in 
a phase that exists only for values ofU/t large enough for a 
Mott insulator to exist at half-filling. This pseudogap is how¬ 
ever very different from the Mott gap. [15, 22, 50] 

AFM We did not systematically study the effect of doping 
on the AFM phase. Nevertheless, Table I is suggestive. Hole¬ 
doping pushes the AFM phase to lower pressures ( higher in¬ 
teraction strengths U/t). For example, at 10% hole-doping, 
the critical pressure for the Neel transition at T/t = 1/20 is 
decreased by about 35% (Fig. 5(c)) compared to the half-filled 


case. By contrast, electron-doping brings the Neel transition 
to higher pressure values (to lower interaction strengths U/t 
) (not shown here). A calculation of the Lindhard function 
shows that this can in part be attributed to better nesting in the 
electron-doped case. 


n 

1.10 

1.00 

0.99 

0.90 

U/t 

3.20 zb 0.05 

4.465 zb 0.005 

4.60 zb 0.05 

6.9 zb 0.1 


TABLE I. Filling dependence of Un/ t for T/t = 1/20. For U/t > 
Pat/ t the AFM phase is stable at that temperature 

The SC State Fig. 5 illustrates the dramatic effects of dop¬ 
ing on the SC phase. For 1% hole-doping (Fig. 5(b)), the 
SC dome is extended by a factor of six on the pressure axis 
at T/t = 1/60 compared with the half-filled case shown in 
Fig. 5(a), while TJ^ (the maximum Tc) is enhanced by ap¬ 
proximately 25%. Suppressing the Mott insulating phase by 
doping allows the SC state to extend its stability far beyond 
the critical interaction strength for the Mott transition at half¬ 
filling (Pm/t) • Comparison of Fig. 5(b) and Fig. 5(c), reveals 
that increasing hole-doping moves to lower pressures. At 
10 % hole doping, decreases slightly (by about 1%) com¬ 
pared to the half-filled case. By contrast. Fig. 5(d) shows that 
electron-doping (10%) displaces to significantly lower 
temperature compared to half-filling (about 15%). The value 
of as a function of doping appears in Table II. As sug¬ 
gested by this Table and Figs. 5(b), 5(c), the maximum dop¬ 
ing for TJ^ is found for intermediate hole-doping at a value 
of U/t that is doping dependent. This is illustrated in Figs. 6 
and 7. Too much hole-doping (10%) or electron-doping (10%) 
effectively reduces TJ^ . Nevertheless, every doping that we 
studied exhibits an enhancement of the range where the SC 
state appears on the pressure axis when compared with half¬ 
filling. 


n 

(ur/t-, PD , t' = o.4i 

(C/r/i; PD , t' = o.8i 

1.10 

(8.3 zb 0.2; 29.0 zb 0.2) 

— 

1.00 

(6.15 zb 0.02; 25.0 zb 0.5) 

(7.78 zb 0.02; 35.0 zb 0.5) 

0.99 

(6.6 zb 0.2; 20.5 zb 0.5) 

— 

0.98 

(6.8 zb 0.2; 18.9 zb 0.1) 

— 

0.97 

(7.2 zb 0.3 ; 19.2 zb 0.5) 

— 

0.96 

(7.3 zb 0.2 ; 19.5 zb 0.5) 

— 

0.94 

(8.0 zb 0.5; 20.8 zb 0.3) 

— 

0.92 

(9.0 zb 0.5; 22.5 zb 0.5) 

— 

0.90 

(10.5 zb 0.5; 25.2 zb 0.2) 

(14.5 zb 0.5; 37.0 zb 0.5) 


TABLE 11. Estimated CDMET values of the maximum superconduct¬ 
ing transition temperature and corresponding interaction strength 
fU^/t\ (TJ^ = 1//^J^) for various fillings and two values 
of frustration t'/t. 

The general aspect of the dome near the triangular right 
part of Figs. 5(a) and 5(b) shows that this section of the SC 
state for the 1% hole-doped case is continuously connected 
to the SC phase at half-filling. For the range of U/t that 
is insulating at half-filling, the SC transition temperature as 
a function of doping at fixed U seems however to vanish 
extremely steeply between 1% and 0% doping, by analogy 
with cuprates. [33, 39, 63] Larger cluster calculations for the 
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0.90 0.95 1.00 1.05 1.10 

n (filling ) 


FIG. 6. Superconducting as a function of filling and U/t pro¬ 
jected in the T — n plane. The actual values of T, f/, n are listed in 
Table II. 



FIG. 7. Superconducting phase diagram combining 'mT — U — n 
space the constant-doping T — U planes of Fig. 5. The line of 
in Fig. 6 also appears on this plot. 


cuprates [78] suggest the existence of a smooth maximum. 
The overall qualitative shape of the SC phase diagram as a 
function oiT — U — n shown in Fig. 7 combines the results of 
Fig. 5. The line of from Fig. 6 also appears on this plot. 
The constant doping scans of this plot are complementary to 
the constant U/t scans of Refs. 35 and 36. 

The order parameter and the uniform susceptibility in the 
SC state (discussed below) illustrate in more detail the evolu¬ 
tion from half-filling to finite doping. 

Order parameter and magnetic susceptibility in the SC state 
Fig. 8 displays the d-wave SC order parameter {dsc ) calcu¬ 
lated using Eq. 10 for = 0.4t. As seen in Fig. 8(a) for 
half-filling, dsc is largest near the first-order transition to 
the insulating phase. This confirms previous T = 0 results 
[74]. Also, as one would expect, as the temperature is raised, 
the magnitude of dsc and the range of pressure where it is 
non-zero decrease. Furthermore, dsc obtained at T = 0 in 
Ref. [74] and at low T in Fig. 8(a) follows qualitatively the 
same pressure dependence as Tc obtained here. 

Fig. 8(b) indicates that the pressure dependence of dsc 
at low temperatures also follows qualitatively that of Tc for 
all other fillings studied, except for the 1% hole-doped case, 
where an anomaly is present near the pseudogap to metal tran¬ 


sition (Fig. 8(c)). This is analog with the unfrustrated square 
lattice where the doping dependence of the low T value of 
dsc does not follow that of Tc in the under-doped regime 
[23, 39]. 

The pseudogap to metal transition leaves some traces in the 
SC state at 1% doping through signatures in certain observ¬ 
ables, such as the uniform susceptibility, = = = 

Eq. 6. Indeed, as shown in Fig. 8(c), near the value 
U = Ummt where the transition between pseudogap and 
metal phases occurs in the normal state, the uniform suscepti¬ 
bility in the SC phase (blue squares) shows a large variation, 
suggesting a crossover in the fundamental metallic proper¬ 
ties of the system, even in the SC state: over a small range 
of U/t (0.3), varies by a factor of 4.3 when Ummt is 
crossed. Upon decreasing pressure, dsc (green diamonds) 
also starts to weaken. This last characteristic is reminiscent of 
the half-filled case, where dsc dies out at the first-order tran¬ 
sition to the insulating phase (Fig. 8(a)). Note by comparing 
Fig. 5(b) and Fig. 8(c) that Tc does not change as drastically 
as Xz{q = O.uj = 0) or dsc upon crossing U = Ummt • The 
value U^, corresponding to the maximum Tc, differs slightly 
from the value where dsc is maximum. 

For the 10% hole-doped case (and 10% electron-doped 
case, not shown), dsc is anti-correlated to . This is illus¬ 
trated in Fig. 8(d). This behavior is expected since it becomes 
more difficult to align the spin of the electrons along a mag¬ 
netic field when more singlet pairs are formed. By contrast 
with the 1% doping case, for 10%, the first-order pseudogap 
to metal transition is not found at accessible T and U/t , so 
X;z shows no peculiar behavior. 

The T = 0 extrapolations of our results are consis¬ 
tent with the variational Monte Carlo results of Watanabe et 
al [30] for the same model. In the normal state they find a 
rapid crossover between two different metallic states at finite- 
doping. They also find that the SC phase for t' = 0.4t is stable 
between U/t = 5 and U/t = 30 for a hole-doping of 8.3%. Our 
numbers dLVoU/t = 6.70 and U/t = 30 for 10% hole-doping at 
T/t = 1/60. 


B. t' = 0.8t 

We move to the case t' = 0.8t, contrasting the results 
with the less frustrated case t' = 0.4t just considered. We 
expect that for larger frustration, both SC and AFM will be 
negatively affected, leaving more room to effects related to 
the Mott transition. However, as we shall see, the effects of 
frustration are much stronger on the AFM than on the SC and 
normal states. 

For compounds with t' = 0.4f, Fig. 9(a) shows that T = 0 
AFM order would mask the Mott transition, and leave either a 
small region of SC at large pressure and very low temperature, 
or coexisting AFM and SC phases. 

Figs. 10(a) and 9(b) at half-filling demonstrate that larger 
frustration, t' = 0.8t, is more detrimental to AFM order than 
to SC order. The maximum Neel transition temperature for 
the AFM phase occurs at (T/t ^ 1/3.5, U/t 8.2) for inter¬ 
mediate frustration, t' = OAt Fig. 9(a), whereas for t' = O.St 








10 


u/t u/t 




(a) (b) 


U/t 



U/t 



(c) 


(d) 


FIG. 8. (color online) Order parameter (dsc ) and uniform susceptibility (xz ) for the superconducting phase of the Hubbard model on the 
anisotropic triangular lattice with t' = 0.4t. Each filling has its specific color and plotmarker in the various panels of this figure and those of 
Fig. 11. Lines are guides to the eye. 

(a) dsc for the half-filled case. Two different temperatures are plotted. 

(b) dsc for various fillings. The temperature is T/t = 1/60. The dashed purple line is an extrapolation. 

(c) dsc (green diamonds) and Xz (blue squares) for the 1% hole-doped case, x^ is divided by 100 to fit on the same vertical axis. Ummt 
denotes the vicinity of interaction strengths where the pseudogap to metal first-order transition discussed in Sec. IV is found. stands for 
the critical interaction strength associated with , the maximum Tc. 

(d) dsc (orange triangles) and Xz (blue squares) for the 10% hole-doped case, x^ is divided by 100 to fit on the same vertical axis. has 
the same significance as in (c). The orange curve for dsc is the same as in (b). 


in Fig. 9(b) one finds (T/t - 1/8.25, U/t - 10). Thus, 
the maximum AFM transition temperature is decreased by a 
factor of about 2.4. By contrast, the maximum Tc decreases 
by about only 25 %. This is the expected effect of frustra¬ 
tion and it agrees qualitatively with FLEX calculations [79] . 
However, FLEX is not valid across a Mott transition. Con¬ 
cerning the nature of the AFM phase, T = 0 studies with 
CDMFT [74] on 2 X 2 cluster, find that between t'/t = 0.7 
and t'/t = 0.8, superconductivity at the first-order transition 
changes from coexisting with a commensurate AFM phase to 
coexisting with a phase that is not magnetically ordered. Vari¬ 


ational Monte Carlo studies find that commensurate AFM is 
stabilized for [76] t'/t < 0.75 or even t'/t < 0.9 [75], con¬ 
sistent with our results. The latter early study however finds 
that magnetic states are always more stable than superconduc¬ 
tivity. 

Long-wavelength AFM fluctuations are also detrimental 
to long-range order since the Mermin-Wagner-Hohenberg the¬ 
orem requires that the Neel transition temperature (T/v) van¬ 
ish in the absence of coupling to the third dimension. The 
T/v lines that we find here are only indicators for the on¬ 
set of the renormalized classical regime where low-frequency 
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U/t 


35. 15. 9.2 6.7 5.3 4.4 



U/t 
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(a) 


(a) 


U/t 


35. 15. 9.2 6.7 5.3 4.4 



(b) 


FIG. 9. Phase diagrams for the Hubbard model on the anisotropic 
triangular lattice. 13 is the inverse temperature in 1/t units. The 
value t — 0.044 eV is used to convert to physical units. [64] Lines 
are guides to the eye. 

(a) Phase diagram at half-filling for an anisotropy parameter t' — 
0.4/. This figure is the same as Fig. 5(a) except that it shows the 
complete antiferromagnetic phase. 

(b) Phase diagram at half-filling for an anisotropy parameter t' — 
0.8/. This figure is the same as Fig. 10(a) except that it shows the 
complete magnetic phase. 


long-wavelength AFM fluctuations become important. Fur¬ 
thermore, as discussed in V A, hole-doping also suppresses 

Tn. 

The above considerations suggest that for /' = 0.8/ our re¬ 
sults for the normal state pseudogap to metal transition and 
for the SC phase are observable at finite temperature in real 
materials. The crossover discussed in Secs. V A and IV also 
occurs for /' = 0.8/ and is thus a fundamental feature of sys¬ 
tems near half-filling. For 1% hole-doping, we found the first- 
order transition at /7// ^ 8. However, the pseudogap to metal 


U/t 

31. 19. 14. 11. 8.9 



(b) 


FIG. 10. (color online) Phase diagrams for the Hubbard model on 
the anisotropic triangular lattice with anisotropy parameter /' = 0.8/ 
and two fillings. The antiferromagnetic state was studied at half¬ 
filling and also forn = 0.9. In the latter case, no AFM was found 
for T down to 1/60 and S < U/t < 30. The value / = 0.044 eV is 
used to convert to physical units. [64] Lines are guides to the eye. 

(a) Phase diagram for the half-filled case. Metallic and insulating 
phases coexist in the blue region. A first-order jump in the normal- 
state double occupancy is found at the blue points. The value of the 
critical point for the Mott transition is (U/t = 7.93, (3 = 10.00). 
With increasing pressure, we find at the red triangles that the Neel 
AFM order parameter disappears (to the right of the red line) through 
a first-order jump for /3 > 10 and through a second-order transition 
for 13 < 10. When pressure is decreased, the paramagnetic metal is 
unstable to the AFM insulator along the green line through a first- 
order jump (^ > 10) or through a second order transition {{3 < 10). 
The orange region is where superconductivity manifests itself. The 
SC state gives way to the insulating phase along a first-order jump 
upon decreasing pressure, and to the metallic phase along a second- 
order transition line upon increasing pressure. 

(b) Phase diagram for the 10% hole-doped case. The transition from 
the SC state to the normal state is second-order. 
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transition and especially the SC phase were particularly diffi¬ 
cult to study extensively due to a worse sign problem, hence 
they are not displayed. 

Results for the 10% doped case are shown in Fig. 10(b). 
The SC phase is present for a much broader range of pres¬ 
sure than in the half-filled case (Fig. 10(a)). Furthermore, the 
SC dome is shifted to lower pressure. These effects of hole¬ 
doping for t' = O.St are very similar to those for t' = OAt. 
Again, the effect of frustration on the superconducting Tc is 
much smaller than on Tat. Indeed, Tn changes from a finite 
value at t'/t = 0.4 (Fig. 5) to zero at t'/t = 0.8, while the 
maximum Tc decreases only by about 30%. 

SC order parameter Fig. 11 displays dsc calculated with 
Eq. 10 for half-filling and for 10% doping. The qualitative 
observations made for t' = 0.4t for 10% doping still hold here 
at larger frustration, namely the pressure dependence of dsc 
at low temperatures follows qualitatively that of Tc and doping 
increases drastically the range of the SC dome on the pressure 
axis. Here too, the uniform susceptibility is anti-correlated to 
dsc ■ 


Additional comparisons with earlier results 

The range of U/t where the SC phase appears for 10% 
hole-doping and t' = O.St, namely 10 < f//t < 26 for 
T/t = 1/60, is similar to that found [30] with the varia¬ 
tional Quantum Monte Carlo method for 8.3% hole-doping 
at T = 0, namely 10 < f//t < 25. 

At half-filling, the shape of the coexistence region in our 
normal state phase diagram differs from that found on a three 
site cluster with CDMFT with exact diagonalization solver by 
Liebsch et al. [64] Also, on the doped three-site cluster the 
pseudogap is not observed in the isotropic limit [80]. These 
differences are not surprising since the ground state entropy of 
clusters with an odd or even number of sites is very different. 
The four site cluster results of Ohashi and al [65] for the 
transition line have a slope of the same sign as us in the T — (7 
plane but the actual values of U and especially T differ, a 
difference that may come from the systematic imaginary-time 
discretization errors of the Hirsch-Fye algorithm. 

Our results extrapolated to T = 0 are consistent with 
earlier CDMFT results obtained with an exact-diagonalization 
impurity solver. [74] Here we also find with decreasing pres¬ 
sure that the system changes from a paramagnetic metal, to 
a superconductor to an AFM insulator. In Ref. 74, the AFM 
phase boundary at sufficiently large frustration coincides with 
the Mott transition. This is consistent with the extrapolation to 
T = 0 of the coincidence in Fig. 10(a) between two transitions 
found with increasing pressure, namely the first-order jump 
between the AFM state and the paramagnetic metal and (when 
AFM long-range order is forbidden) the first-order jump be¬ 
tween the Mott insulator and the paramagnetic metal. At 
smaller frustration, this does not occur: Upon decreasing pres¬ 
sure (or increasing U), the AFM transition occurs before the 
Mott transition, as can be seen for t'/t = 0.4 in Fig. 9. Also, 
as seen in Fig. 9, the critical value of U/t for the Mott tran¬ 
sition at half-filling increases with frustration t'/t, in agree¬ 


ment with Refs. [74, 81]. This is reflected in the qualita¬ 
tive shape of the magenta region in Fig. I. The Mott transi¬ 
tion on the anisotropic triangular lattice has been studied with 
many other methods, for example path integral renormaliza¬ 
tion group [82]. The critical value of U/t found with the lat¬ 
ter method is smaller than that in CDMFT. Note finally that 
it has been found earlier in different contexts that depending 
on frustration t'/t, AFM transitions can be of first or second 
order [83-85]. 

The T = 0 magnetic phases found for various values of 
t' /t and U/t with the variational cluster approximation [86- 
88 ] exact diagonalization of an effective model [89], varia¬ 
tional Monte Carlo, [75, 76] and dual-fermions [87] are sim¬ 
ilar to those found with CDMFT in Ref. 74, except that here 
and in Ref. 74 the possibility of spiral order has not been 
investigated. Nevertheless, according to Refs. [75, 76, 90], 
commensurate AFM of the type we find can be stable up to 
t'/t 0.8. In [75] superconductivity is always less stable 
than AFM but this is not so in weak correlation approaches 
such as FLEX [79] or functional renormalization group [91]. 
With the latter approach, superconductivity is also obtained 
for t' /t = 1 but symmetry considerations in this case are dif¬ 
ferent. [92] At half-filling, it has been widely appreciated for 
a long time by many methods that superconductivity is sta¬ 
ble for a wide range of values of t'/t\ in variational [93, 94] 
CDMFT [74] Gutzwilier [95] resonating valence bond [96] 
approaches. 


VI. DISCUSSION 
A. Broken-symmetry phases 

Although phase transitions in CDMFT are renormal¬ 
ized by local dynamical fiuctuations, they essentially have a 
mean-field character. In particular, they do not satisfy the 
Mermin-Wagner-Hohenberg theorem that forbids continuous- 
symmetry breaking at finite temperature in two dimen¬ 
sions. [97] This is especially relevant for AFM order, which 
we overestimate. On the other hand, the superconducting 
critical temperature Tc found in our phase diagrams physi¬ 
cally represents the dynamical mean-field transition temper¬ 
ature T^ below which Cooper pairing occurs locally in the 
cluster. [39] Long-wavelength quantum and thermal fiuctua¬ 
tions in the amplitude and phase of the order parameter dsc 
should lead to an actual Kosterlitz-Thouless transition temper¬ 
ature smaller than T^ [98-100] With increasing cluster sizes, 
the dynamical mean-field Tc have been shown to converge to 
a finite value on the square lattice. [27] 


B. Strongly correlated superconductivity, superconducting 
dome and AFM quantum critical point 

The link between antiferromagnetic quantum critical point 
and unconventional superconductivity is well documented, es¬ 
pecially in the field of heavy-fermion materials. [101] Nu¬ 
merical simulations with methods very close to those used 
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FIG. 11. (color online) The order parameter (dsc ) and uniform sus¬ 
ceptibility (xz ) in the superconducting phase of the Hubbard model 
on the anisotropic triangular lattice with frustration value t' = 0.8t. 
Each doping has the same specific color and plotmarker as in Fig. 8. 
Lines are guide to the eye. 

(a) dsc for the half-filled case. Two different temperatures are plot¬ 
ted. 

(b) dsc for various fillings and Xz for the 10% hole-doped case. The 
temperature is T/t = 1/60. stands for the superconducting crit¬ 
ical interaction strength associated with . The spin susceptibility 
Xz calculated with Eq. 6 in dimensionless units is divided by 200 to 
fit on the same vertical axis. 


here confirm this intimate connection for the Anderson lat¬ 
tice model of heavy-fermions: [102] Indeed, one finds a su¬ 
perconducting dome that systematically surrounds the antifer¬ 
romagnetic quantum critical point. The same type of con¬ 
nection to a quantum-critical point has been proposed for 
cuprates. [103, 104] We suggest that this connection between 
AFM quantum critical point and maximum Tc is present when 


the interaction strength is not large enough to lead to a Mott 
transition. In that case pairing occurs through the exchange of 
long wavelength antiferromagnetic fluctuations. [105-108] 

The top panel of Fig. 9 shows that for half-fllled organics, 
where a Mott transition is clearly observed, superconductiv¬ 
ity is near the Mott transition, not near the antiferromagnetic 
quantum critical point. And as we dope. Fig. 5(b) shows that 
the superconducting dome surrounds the pseudogap to metal 
transition that is the flnite-doping remnant of the Mott transi¬ 
tion. 

A schematic phase diagram is displayed in Fig. 12. In the 
normal state, there is a flrst-order phase transition whose co¬ 
existence region is represented in blue. The maximum of 
the superconducting transition temperature is controlled by 
the opening of the pseudogap, namely by the position of the 
flrst-order phase transition or its continuation, not by the an¬ 
tiferromagnetic quantum critical lines at the end of the anti¬ 
ferromagnetic three-dimensional dome delimited by the red 
region: even in the absence of the antiferromagnetic phase, 
superconductivity survives, as can be verifled from Fig. 10b. 
This is a characteristic of strongly-correlated superconductiv¬ 
ity in doped Mott insulators. We stress however that at large 
doping we observe crossovers in the normal state, but we can¬ 
not calculate at low enough temperature to confirm if the flrst- 
order transition survives. It can in principle be replaced by a 
T = 0 second order transition line or disappear at a critical 
point. 

The mechanism for superconductivity in the organics is 
thus clearly different from that associated with an antifer¬ 
romagnetic quantum critical point. This type of strongly- 
correlated superconductivity is controlled by short-range 
AFM correlations, namely superexchange J = jU, as 
found early on in slave-boson calculations [109] and more re¬ 
cently with CDMFT [61]. This is reviewed in Ref. 15. 

Further studies will be needed to clarify the detailed cause 
of the Tc dome. Consider the optimal Tc in Figs. 10(b), 5(b), 
5(c), 5(d). The decrease of Tc from optimal towards low pres¬ 
sure (on the left) scales like J (i.e. like a straight line). On 
the other hand, the opening of a pseudogap could also explain 
this decrease because a pseudogap in the density of states at 
the Fermi level leaves fewer states to pair. 


C. Contact with experiment 

Our phase diagrams at half-filling (Figs. 5(a), 10(a)), re¬ 
veal interesting similarities with the experimental phase di¬ 
agrams of the BEDT compounds but even more so with 
the Pd(dmit) 2 compounds. [3, 4] In the temperature-pressure 
plane, the shape of the region where the SC phase exists and 
the shape of the metal-insulator coexistence region are in qual¬ 
itative agreement with experiment. The maximum SC transi¬ 
tion temperature T^ as a function of pressure coincides with 
the phase transition to the insulator, also like in experiment. 

We And that as frustration is increased, both the max¬ 
imum T/v and T^ decrease, the decrease in the tendency to 
magnetic order being however much more pronounced. These 
observations are consistent with experiments where com- 
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FIG. 12. Schematic phase diagram for layered organics. Orange 
planes indicate the three-dimensional region where superconductiv¬ 
ity is present. The antiferromagnetic three-dimensional dome is be¬ 
low the red lines. In the normal state, there is a surface of first- 
order transition whose coexistence region is delimited by the blue 
lines. This first-order transition separates a pseudogap phase with 
small double-occupancy from a correlated metal with larger double¬ 
occupancy, by analogy with Ref. 35. At larger U and doping, the 
sign problem prevents us from observing directly the first-order tran¬ 
sition or its continuation as a different phenomenon. Frustration can 
appreciably displace and eventually completely eliminate the antifer¬ 
romagnetic region. [74, 110] 


pounds with higher frustration t'/t have lower Tc compared 
to less frustrated materials. For example, /^-(ET) 2 Cu 2 (CN )3 
(t' = O.St) and /^-(ET) 2 Cu[N(CN) 2 ]Cl (t' = OAt) dis¬ 
play of respectively 3.9K and 13. IK. Also, as frustration 
is increased in real compounds, AEM order is greatly sup¬ 
pressed. Eor instance, /^-(ET) 2 Cu[N(CN) 2 ]Cl has AEM order 
at low temperatures and pressures, while /^-(ET) 2 Cu 2 (CN )3 is 
a spin-liquid candidate. [1,7] 

Eig. 10 at = O.St accounts for the experimen¬ 
tal results on the SC phase both at half-filling [2] and at 
10% doping. [16] Indeed, in half-filled /^-(ET) 2 Cu 2 (CN )3 a 
superconductor-insulator phase transition is observed upon 
decreasing pressure, as in Eig. 10(a). On the other hand 
/^-(ET) 4 Hg 2 . 89 Br 8 , the 10% hole-doped counterpart of k- 
(ET) 2 Cu 2 (CN )3 , presents a dome shaped [16, 20] Tc similar 
to Eig. 10(b) except that the dome is more asymmetrical in 
experiment. 

Eig. 10 for t' = O.St also shows that the range of pres¬ 
sure where SC appears at 10% hole-doping and T/t = 1/60 
is multiplied by about nine compared to the half-filled case. 
Rough extrapolation of the superconducting dome to T = 0 
gives a range of pressure that increases by a factor four to six 
in going from half-filling to 10% doping. Experimentally the 
dome is extended by approximately six for the same value of 
frustration. Indeed, in half-filled /^-(ET) 2 Cu 2 (CN )3 , super¬ 
conductivity occurs over a range of 0.25 GPA according to 
Ref. [Ill] while Ref. [16] finds superconductivity for a range 
of 1.5 GPa in the doped compound /i:-(ET) 4 Hg 2 . 89 Br 8 . 

TJ^ for /^-(ET) 2 Cu 2 (CN )3 is 3.9i^ while it is about 7K for 
the doped counterpart, a factor of 1.8. Our results show a very 
slight decrease of at 10% hole-doping. However, is 
increased for intermediate doping (Table II, Pig. 6) compared 


to half-filling. While that increase for intermediate doping is 
proven for an anisotropy parameter t' = 0.4t, the case t' = 
0.8/ should be similar. 

Oike et al. [16] and Taniguchi et al. [20] also found in 
the Hall coefficient of /^-(ET) 4 Hg 2 . 89 Br 8 a rapid crossover 
around 0.5 GPa. This pressure corresponds closely to the 
maximum of the superconducting transition temperature [16] 
. Although for /'// = 0.8 and n = 0.9 we cannot reach 
low-enough temperature to detect the first-order pseudogap to 
metal transition, the /'// = 0.4 low-doping results of Pig. 5(b) 
strongly suggest that is controlled by that transition. On 
the square lattice, one finds analogous results: [39] A low tem¬ 
perature first-order transition between a pseudogap metallic 
phase and a strongly-correlated metal ends at a critical point 
above which a line of crossovers appears. [22, 23] This line 
of crossovers is a Widom line, a general phenomenon found 
in the supercritical region of first-order transitions. [66] Re¬ 
markably, is near the intersection of the superconducting 
dome and of the Widom line. [21, 39] This leads us to impor¬ 
tant predictions for experiment. 


D. Predictions 

The first-order pseudogap to metal transition, observed 
theoretically on the square lattice, [22, 35, 36] has also been 
seen as a sharp crossover in larger cluster calculations in the 
Dynamical Cluster Approximation [112] and in variational 
Quantum Monte Carlo. [113] Our work shows that this tran¬ 
sition also occurs on the anisotropic triangular lattice. It had 
been observed before as a sharp crossover. [30] The experi¬ 
mental results on the doped BEDT /^-(ET) 4 Hg 2 . 89 Br 8 [16, 20] 
can thus be interpreted as an observation of the crossovers as¬ 
sociated to this pseudogap to metal transition. Based on our 
phase diagrams, we predict that in yet to be synthesized very 
low-doping organic materials with pressure-induced transi¬ 
tions, remnants of this transition could be detectable in the 
SC state. It would manifest itself via observables such as the 
uniform susceptibility, the SC order parameter or by a strong 
crossover of many properties as a function of pressure near 

rp m 
c 

We also predict that in electron-doped compounds, Tc is 
decreased but the range of SC on the pressure axis is still in¬ 
creased compared to the half-filled case. 

Einally, our results also suggest that for frustration high 
enough that magnetically ordered phases are absent, the nor¬ 
mal state underlying the SC state of the doped compounds 
should display a first-order transition between a pseudogap 
and a more metallic state at sufficiently low temperatures. 
This is our most important prediction. As usual the normal 
state can be revealed by applying a magnetic field. Experi¬ 
ments at half-filling [114] on /i:-(BEDT-TTE) 2 Cu[N(CN) 2 ]Cl 
suggest that the magnetic fields necessary to destroy the SC 
state are easily accessible. 
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E. Limitations and perspectives 

The broad agreement that we find with experiment and, in 
the low temperature limit, with variational Monte Carlo meth¬ 
ods, [30] suggest that the important physics in the layered 
organics arises from strong on-site repulsion U and nearest- 
neighbor superexchange J. Increasing the cluster size would 
allow a better representation of the long-wavelength fiuctua- 
tions beyond mean-field theory. It would produce more quan¬ 
titative phase boundaries, as discussed in Sec. VIA above. 
But since the existence of the superconducting phase itself has 
already been established by finite-size studies on the square 
lattice, [27] it is highly unlikely that larger cluster studies 
would change this. 

The continuous time Quantum Monte Carlo impurity solver 
in the hybridization expansion that we have used here (CT- 
HYB) is for now the only Monte Carlo approach that allows 
calculations in the range of large U/t and frustration needed 
for the layered organics. Recall that U/t = 14 at the maxi¬ 
mum of the superconducting dome for f'/f = 0.8 and that U/t 
is as large as 30 at the lowest temperature end of the dome. 
Expansions in power ofU /t that are used as Quantum Monte 
Carlo impurity solvers for larger clusters [45] fail for such 
large values of U/t and frustration because of a sign prob¬ 
lem and because of the large expansion orders that are needed. 
Even on the 4-site cluster that we use, the average sign is im¬ 
measurably small with Rubtsov’s algorithm fox U/t = 14.5 
and T/t = 1/40, namely near the maximum of the super¬ 
conducting dome for t'/t = 0.8. [45] Eor high-temperature, 
(T/t > 0.06) the cluster-size dependence has been shown to 
be negligible. [115] While CT-HYB provides a method to ac¬ 
cess the large values of U/t and t'/t that we need, the compu¬ 
tation time increases exponentially with system size, making 
larger clusters unfeasible with present resources. The same 
size limitation applies to exact-diagonalization solvers [52] 
that, in addition, rely on a finite bath, by contrast with CT- 
HYB where the bath is infinite. A method has recently been 
proposed [116] to increase bath size in exact diagonalization 
solvers but extensions to higher temperatures have not been 
tested and implementations on clusters have not been done 
yet. 

Although AEM and SC were considered on equal footing, 
we neglected the possibility of non-vanishing order parame¬ 
ters for both AEM and SC simultaneously. This might oc¬ 
cur in some regions at finite doping where we found that both 
AEM and SC separately can develop long-range order. The 
question of simultaneous non-vanishing order parameters is 
an interesting question but it is a detail at this early stage of 
investigations. The fact that for t'/t = 0.8 model of the doped 
organic /^-(ET) 4 Hg 2 . 89 Br 8 there is a superconducting dome 
far from the AEM phase is one of the crucial proofs that the 
maximum of the dome does not come from an AEM quantum 
critical point. 

So far, studies of the Hubbard model have shown that they 
are capable of capturing essential features of materials such as 
cuprates, /t:-(ET) 2 X or Pd(dmit) 2 , even if that model neglects 
some physical effects. Additional calculations taking into ac¬ 
count Coulomb interaction between nearest neighbors (with 


the extended Hubbard model [38, 117]), electron-phonon in¬ 
teractions [118], or the third spatial dimension [119, 120] 
would allow one to capture increasing details of these fasci¬ 
nating compounds, but the overall agreement with experiment 
that has been found suggests that the local and near-neighbor 
superexchange aspects of the Hubbard model capture the es¬ 
sential physics. 

Eurther investigations of the normal state pseudogap and 
of properties of the strongly correlated superconducting phase 
are planned. 


VII. CONCLUSION 

Based on CDMET calculations of normal, superconduct¬ 
ing and antiferromagnetic phase diagrams for the Hubbard 
model on the anisotropic triangular lattice, we arrive at the 
following conclusions. These phase diagrams are very sim¬ 
ilar to experimental observations, both at half-filling and at 
finite doping. Upon doping, superconductivity is enhanced; 
in particular it occurs over a much broader range of pressures 
(U/t). Smaller anisotropy, or equivalently larger frustration 
(t'/t ^ 1), diminishes antiferromagnetic and superconducting 
transition temperatures but antiferromagnetism is much more 
affected. 

In the normal state, a first-order pseudogap to metal transi¬ 
tion occurs at finite doping and low temperature. The transi¬ 
tion is continuously connected to the Mott transition at half¬ 
filling, as on the square lattice case, [35] yet the pseudogap 
is different from the Mott gap. The decrease in the density 
of states at the Eermi level associated with the pseudogap and 
the decrease of superexchange J when pressure decreases can 
both contribute to the decrease of Tc with decreasing pressure. 
We claim that competing antiferromagnetism is not an expla¬ 
nation for the dome in the doped organics. 

We predict that for very lightly hole-doped compounds, the 
pseudogap to metal transition leaves some subtle traces in the 
superconducting state. Our most important prediction is that 
the normal state that will be revealed by destroying the super¬ 
conducting state with a magnetic field in lightly-doped highly- 
frustrated compounds will show this first-order transition be¬ 
tween two conducting phases, one with a pseudogap and the 
other one metallic. It is this transition that should control the 
crossovers at finite temperature as well as the location of the 
maximal superconducting transition temperature. Einally, we 
predict that electron-doping should lead to a reduced maxi¬ 
mum Tc. 

Eor a model of heavy fermions solved with the same set of 
methods as those used here, the maximum of the supercon¬ 
ducting dome can be correlated with the location of an anti¬ 
ferromagnetic quantum critical point for interaction strengths 
that are not large enough to lead to a Mott transition. [102] 
In our case, the superconducting dome in the doped organ¬ 
ics surrounds the finite-doping extension of the zero-doping 
first-order Mott transition (whenever it is directly observable). 
It is definitely not attached to an antiferromagnetic quantum 
critical point, a conclusion that can also be verified exper¬ 
imentally. This result should clearly impact understanding 
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of strongly-correlated superconductivity in all doped Mott in¬ 
sulators; not only layered organic superconductors, but also 
high-temperature cuprate superconductors. 

ACKNOWLEDGMENTS 


nais, and R. Nourafkan for suggestions and detailed comments 
on the manuscript. This work was supported by the Nat¬ 
ural Sciences and Engineering Research Council of Canada 
(NSERC) and the Tier I Canada Research Chair Program (A.- 
M.S.T.). Simulations were performed on computers provided 
by CEI, MELS, Calcul Quebec and Compute Canada. 


We are grateful to G. Sordi for numerous discussions, to 
Wei Wu for some calculations and to M. Gingras, C. Bourbon- 


[1] S. Lefebvre, P. Wzietek, S. Brown, C. Bourbonnais, 

D. Jerome, C. Meziere, Cziere, M. Fourmigue, and P. Batail, 
Physical review letters 85, 5420 (2000). 

[2] Y. Shimizu, H. Kasahara, T. Furuta, K. Miyagawa, K. Kanoda, 
M. Maesato, and G. Saito, Phys. Rev. B 81, 224508 (2010). 

[3] Y. Shimizu, H. Akimoto, H. Tsujii, A. Tajima, and R. Kato, 
Phys. Rev. Lett. 99, 256403 (2007). 

[4] J.-I. Yamaura, A. Nakao, and R. Kato, Journal 
of the Physical Society of Japan 73, 976 (2004), 
http://dx.doi.org/10.1143/JPSJ.73.976. 

[5] P. Limelette, P. Wzietek, S. Florens, A. Georges, T. A. Costi, 
C. Pasquier, D. Jerome, C. Meziere, and P. Batail, Physical 
review letters 91, 016401 (2003). 

[6] J. Merino and O. Gunnarsson, Phys. Rev. B 89, 245130 (2014). 

[7] Y Shimizu, K. Miyagawa, K. Kanoda, M. Maesato, and 
G. Saito, Phys. Rev. Lett. 91, 107001 (2003). 

[8] M. Yamashita, N. Nakata, Y Senshu, M. Nagata, H. M. Ya¬ 
mamoto, R. Kato, T. Shibauchi, and Y Matsuda, Science 328, 
1246 (2010). 

[9] B. J. Powell and R. H. McKenzie, Journal of Physics: Con¬ 
densed Matter 18, R827 (2006). 

[10] B. J. Powell and R. H. McKenzie, Reports on Progress in 
Physics 74, 056501 (2011). 

[11] H. Kino and H. Fukuyama, Journal of the Physical Society of 
Japan 65, 2158 (1996). 

[12] R. T. Clay, H. Li, and S. Mazumdar, Phys. Rev. Lett. 101, 
166403 (2008). 

[13] S. Dayal, R. T. Clay, and S. Mazumdar, Phys. Rev. B 85, 
165141 (2012). 

[14] M. Imada, A. Fujimori, and Y Tokura, Rev. Mod. Phys. 70, 
1039 (1998). 

[15] A.-M. S. Tremblay, in Emergent Phenomena in Correlated 
Matter Modeling and Simulation, Vol. 3, edited by E. Pavarini, 

E. Koch, and U. Schollwock (Verlag des Forschungszentrum, 
2013) Chap. 10. 

[16] H. Oike, K. Miyagawa, H. Taniguchi, and K. Kanoda, Phys. 
Rev. Lett. 114, 067002 (2015). 

[17] R. N. Lyubovskaya, E. A. Zhilyaeva, A. V. Zvarykina, V. N. 
Laukhin, R. B. Lyubovskii, and S. 1. Pesotskii, JETP Lett. 45, 
530 (1987). 

[18] R. N. Lyubovskaya, E. I. Zhilyaeva, S. I. Pesotskii, R. B. 
Lyubovskii, L. O. Atovmyan, O. A. D’yachenko, and T. G. 
Takhirov, JETP Lett. 46, 188 (1987). 

[19] S. L. Bud’ko, A. G. Gapotchenko, A. E. Luppov, R. N. 
Lyubovskaya, and R. B. Lyubovskii, Sov. Phys. JETP 74, 983 
(1992). 

[20] H. Taniguchi, T. Okuhata, T. Nagai, K. Satoh, N. Mori, 
Y. Shimizu, M. Hedo, and Y. Uwatoko, Journal 
of the Physical Society of Japan 76, 113709 (2007), 
http://dx.doi.org/10.1143/JPSJ.76.113709. 


[21] O. Cyr-Choiniere, D. LeBoeuf, S. Badoux, S. Dufour- 
Beausejour, D. A. Bonn, W. N. Hardy, R. Liang, N. Doiron- 
Leyraud, and L. Taillefer, ArXiv e-prints (2015), 
arXiv: 1503.02033 [cond-mat.supr-con] . 

[22] G. Sordi, P Semon, K. Haule, and A.-M. S. Tremblay, Scien¬ 
tific Reports 2 (2012), 10.1038/srep00547. 

[23] G. Sordi, P. Semon, K. Haule, and A.-M. S. Tremblay, Phys. 
Rev. B 87, 041101 (2013). 

[24] H. C. Kandpal, I. Opahle, Y.-Z. Zhang, H. O. Jeschke, and 
R. Valenti, Phys. Rev. Lett. 103, 067004 (2009). 

[25] K. Nakamura, Y. Yoshimoto, T. Kosugi, R. Arita, and 
M. Imada, Journal of the Physical Society of Japan 78, 083710 
(2009). 

[26] A. Georges, G. Kotliar, W. Krauth, and M. J. Rozenberg, Rev. 
Mod. Phys. 68 , 13 (1996). 

[27] T. A. Maier, M. Jarrell, T. C. Schulthess, P. R. C. Kent, and 
J. B. White, Phys. Rev. Lett. 95, 237001 (2005). 

[28] G. Kotliar, S. Y Savrasov, K. Haule, V. S. Oudovenko, O. Par- 
collet, and C. A. Marianetti, Reviews of Modern Physics 78, 
865 (2006). 

[29] A. M. S. Tremblay, B. Kyung, and D. Senechal, Low Temp. 
Phys. 32, 424 (2006). 

[30] T. Watanabe, H. Yokoyama, and M. Ogata, JPS conference 
proceedings (2014), 10.7566/JPSCP.3.013004. 

[31] We thank A. Reymbaut for drawing this figure. 

[32] P. Semon and A.-M. S. Tremblay, Phys. Rev. B 85, 201101 

( 2012 ). 

[33] P. Semon, G. Sordi, and A.-M. S. Tremblay, Phys. Rev. B 89, 
165113 (2014). 

[34] P Semon, C.-H. Yee, K. Haule, and A.-M. S. Tremblay, Phys. 
Rev. B 90, 075149 (2014). 

[35] G. Sordi, K. Haule, and A. M. S. Tremblay, Phys. Rev. Lett. 
104, 226402 (2010). 

[36] G. Sordi, K. Haule, and A.-M. S. Tremblay, Physical Review 
B 84, 075161 (2011). 

[37] K. Nakamura, Y. Yoshimoto, T. Kosugi, R. Arita, and 
M. Imada, Journal of the Physical Society of Japan 78, 083710 
(2009), http://dx.doi. 0 rg/lO.l 143/JPSJ.78.083710. 

[38] D. Senechal, A. G. R. Day, V. Bouliane, and A.-M. Tremblay, 
Phys. Rev. B 87, 075123 (2013). 

[39] G. Sordi, P. Semon, K. Haule, and A.-M. S. Tremblay, Phys. 
Rev. Lett. 108, 216401 (2012). 

[40] G. Kotliar, S. Y. Savrasov, G. Palsson, and G. Biroli, Phys. 
Rev. Lett. 87, 186401 (2001). 

[41] T. Maier, M. Jarrell, T. Pruschke, and M. H. Hettler, Reviews 
of Modem Physics 77, 1027 (2005). 

[42] D. Senechal, arXiv preprint arXiv:0806.2690 (2008). 

[43] E. Koch, G. Sangiovanni, and O. Gunnarsson, Phys. Rev. B 
78, 115102 (2008). 

[44] E. Gull, P. Werner, A. Millis, and M. Troyer, Phys. Rev. B 76, 



17 


235123 (2007). 

[45] E. Gull, A. J. Millis, A. L Lichtenstein, A. N. Rubtsov, 
M. Troyer, and R Werner, Reviews of Modern Physics 83, 
349 (2011). 

[46] P. Werner, A. Comanac, L. de’ Medici, M. Troyer, and A. J. 
Millis, Phys. Rev. Lett. 97, 076405 (2006). 

[47] K. Haule, Phys. Rev. B 75, 155113 (2007). 

[48] M. H. Hettler, A. N. Tahvildar-Zadeh, M. Jarrell, T. Pruschke, 
and H. R. Krishnamurthy, Phys. Rev. B 58, R7475 (1998). 

[49] C. Huscroft, M. Jarrell, T. Maier, S. Moukouri, and A. N. 
Tahvildarzadeh, Phys. Rev. Lett. 86 , 139 (2001). 

[50] D. Senechal and A.-M. S. Tremblay, Phys. Rev. Lett. 92, 
126401 (2004). 

[51] B. Kyung, S. S. Kancharla, D. Senechal, A.-M. S. Tremblay, 
M. Civelli, and G. Kotliar, Physical Review B (Condensed 
Matter and Materials Physics) 73, 165114 (2006). 

[52] A. Liebsch and N.-H. Tong, Phys. Rev. B 80, 165126 (2009). 

[53] E. Gull, O. Parcollet, P. Werner, and A. J. Millis, Phys. Rev. 
B 80, 245102 (2009). 

[54] M. Eerrero, P S. Cornaglia, L. De Leo, O. Parcollet, 
G. Kotliar, and A. Georges, Phys. Rev. B 80, 064501 (2009). 

[55] H. Park, K. Haule, and G. Kotliar, Physical Review Letters 
101, 186403 (2008). 

[56] M. Balzer, B. Kyung, D. Senechal, A.-M. S. Tremblay, and 
M. Potthoff, Euro. Phys. Lett. 85, 17002 (2009). 

[57] D. Senechal, P.-L. Lavertu, M.-A. Marois, and A.-M. S. Trem¬ 
blay, Phys. Rev. Lett. 94, 156404 (2005). 

[58] M. Capone and G. Kotliar, Phys. Rev. B 74, 054513 (2006). 

[59] K. Haule and G. Kotliar, Physical Review B (Condensed Mat¬ 
ter and Materials Physics) 76, 104509 (2007). 

[60] M. Aichhorn, E. Arrigoni, M. Potthoff, and W. Hanke, Phys. 
Rev. B 74, 024508 (2006). 

[61] S. S. Kancharla, B. Kyung, D. Senechal, M. Civelli, 
M. Capone, G. Kotliar, and A.-M. S. Tremblay, Phys. Rev. 
B77, 184516(2008). 

[62] E. Gull, O. Parcollet, and A. J. Millis, Phys. Rev. Lett. 110, 
216405 (2013). 

[63] L. Eratino, P. Semon, G. Sordi, and A. M. S. Tremblay, un¬ 
published. 

[64] A. Liebsch, H. Ishida, and J. Merino, Physical Review B 79, 
195108 (2009). 

[65] T. Ohashi, T. Momoi, H. Tsunetsugu, and N. Kawakami, 
Physical review letters 100, 076402 (2008). 

[66] L. Xu, P. Kumar, S. V. Buldyrev, S.-H. Chen, 

P. H. Poole, E. Sciortino, and H. E. Stanley, 
Proc. Nat. Acad. Sci. USA 102, 16558 (2005), 

http://www.pnas.org/content/102/46/16558.full.pdf-i-html. 

[67] G. G. Simeoni, T. Bryk, E. A. Gorelli, M. Krisch, G. Ruocco, 
M. Santoro, and T. Scopigno, Nature Physics 6, 503 (2010). 

[68] H. Alloul, Comptes Rendus Physique 15, 519 (2014). 

[69] The critical line ending the first-order surface plunges to very 
low temperature as U or doping increase. The first-order tran¬ 
sition may end at large U or continue as a quantum-critical 
line. We can only confirm that crossovers exist at high temper¬ 
ature. The sign problem prevents us from exploring very low 
temperature. 

[70] G. Kotliar, S. Murthy, and M. J. Rozenberg, Phys. Rev. Lett. 
89, 046401 (2002). 

[71] A. Macridin, M. Jarrell, and T. Maier, Phys. Rev. B 74, 
085104 (2006). 

[72] E. Khatami, K. Mikelsons, D. Galanakis, A. Macridin, 
J. Moreno, R. T. Scalettar, and M. Jarrell, Phys. Rev. B 81, 
201101 ( 2010 ). 

[73] T. Komatsu, N. Matsukawa, T. Inoue, and G. Saito, Jour¬ 


nal of the Physical Society of Japan 65, 1340 (1996), 
http://dx.doi.org/10.1143/JPSJ.65.1340. 

[74] B. Kyung and A.-M. S. Tremblay, Physical review letters 97, 
046402 (2006). 

[75] T. Watanabe, H. Yokoyama, Y. Tanaka, and J. Inoue, Phys. 
Rev. B 77, 214505 (2008). 

[76] L. E. Tocchio, H. Eeldner, E. Becca, R. Valenti, and C. Gros, 
Phys. Rev. B 87, 035143 (2013). 

[77] On the square lattice, [35, 36] the chemical potential was var¬ 
ied on a fine scale whereas here it is the interaction that is 
varied slowly. The first-order transition occurs on a surface in 
the T — U — /i plane that can be traversed in several directions. 

[78] E. Gull, O. Parcollet, and A. J. Millis, Physical Review Letters 
110, 216405 (2013). 

[79] H. Kino and H. Kontani, Journal of the Physical Society of 
Japan 67, 3691 (1998). 

[80] B. Kyung, Phys. Rev. B 75, 033102 (2007). 

[81] O. Parcollet, G. Biroli, and G. Kotliar, Phys. Rev. Lett. 92, 
226402 (2004). 

[82] H. Morita, S. Watanabe, and M. Imada, Journal 

of the Physical Society of Japan 71, 2109 (2002), 
http://dx.doi.org/10.1143/JPSJ.71.2109. 

[83] Z.-Q. Yu and L. Yin, Phys. Rev. B 81, 195122 (2010). 

[84] H. Kondo and T. Moriya, Journal of the 

Physical Society of Japan 65, 2559 (1996), 

http://dx.doi.org/10.1143/JPSJ.65.2559. 

[85] W. Hofstetter and D. Vollhardt, Annalen der Physik 510, 48 
(1998). 

[86] P. Sahebsara and D. Senechal, Phys. Rev. Lett. 100, 136402 
(2008). 

[87] M. Laubach, R. Thomale, C. Platt, W. Hanke, and G. Li, Phys. 
Rev. B 91, 245125 (2015). 

[88] A. Yamada, Physical Review B 89, 195108 (2014). 

[89] H.-Y. Yang, A. Lauchli, E. Mila, and K. Schmidt, Phys. Rev. 
Lett. 105, 267204 (2010). 

[90] J. Reuther and R. Thomale, Phys. Rev. B 83, 024402 (2011). 

[91] S. W. Tsai and J. B. Marston, Canadian Journal of Physics 79, 
1463 (2001), http://dx.doi.org/10.1139/p01-085. 

[92] C. Honerkamp, Phys. Rev. B 68, 104510 (2003). 

[93] J. Liu, J. Schmalian, and N. Trivedi, Phys. Rev. Lett. 94, 
127003 (2005). 

[94] T. Watanabe, H. Yokoyama, Y. Tanaka, and J. ichiro Inoue, 
Journal of the Physical Society of Japan 75, 074707 (2006), 
http://dx.doi.org/10.1143/JPSJ.75.074707. 

[95] J. Y. Gan, Y. Chen, Z. B. Su, and E. C. Zhang, Phys. Rev. Lett. 
94, 067005 (2005). 

[96] B. J. Powell and R. H. McKenzie, Phys. Rev. Lett. 94, 047004 
(2005). 

[97] N. D. Mermin and H. Wagner, Phys. Rev. Lett. 17, 1133 
(1966). 

[98] 1. Ussishkin, S. L. Sondhi, and D. A. Huse, Phys. Rev. Lett. 
89, 287001 (2002). 

[99] D. Podolsky, S. Raghu, and A. Vishwanath, Phys. Rev. Lett. 
99, 117004 (2007). 

[100] Z. Tesanovic, Nat Phys 4, 408 (2008). 

[101] Y.-f. Yang and D. Pines, Proceedings of the Na¬ 
tional Academy of Sciences 111, 8398 (2014), 

http://www.pnas. 0 rg/c 0 ntent/l 11/23/8398.full.pdf-Hhtml. 

[102] W. Wu and A.-M.-S. Tremblay, Phys. Rev. X 5, 011019 
(2015). 

[103] B. J. Ramshaw, S. E. Sebastian, R. D. McDonald, J. Day, 
B. Tam, Z. Zhu, J. B. Betts, R. Liang, D. A. Bonn, 
W. N. Hardy, and N. Harrison, ArXiv e-prints (2014), 
arXiv: 1409.3990 [cond-mat.supr-con]. 


18 


[104] L. Taillefer, Annual Review of Condensed Matter Physics 
1, 51 (2010), http://dx.doi.org/10.1146/annurev-conmatphys- 
070909-104117. 

[105] M. T. Beal-Monod, C. Bourbonnais, and V. J. Emery, Phys. 
Rev. B 34, 7716 (1986). 

[106] D. J. Scalapino, J. Loh, E., and J. E. Hirsch, Phys. Rev. B 34, 
8190 (1986). 

[107] N. E. Bickers, D. J. Scalapino, and S. R. White, Phys. Rev. 
Lett. 62, 961 (1989). 

[108] T. A. Maier, A. Macridin, M. Jarrell, and D. J. Scalapino, 
Phys. Rev. B 76, 144516 (2007). 

[109] G. Kotliar and J. Liu, Phys. Rev. Lett. 61, 1784 (1988). 

[110] P. Sahebsara and D. Senechal, Phys. Rev. Lett. 97, 257004 
(2006). 

[111] Y. Kurosaki, Y. Shimizu, K. Miyagawa, K. Kanoda, and 
G. Saito, Phys. Rev. Lett. 95, 177001 (2005). 

[112] P. Werner, E. Gull, O. Parcollet, and A. J. Millis, Phys. Rev. 


B 80, 045120 (2009). 

[113] H. Yokoyama, M. Ogata, Y. Tanaka, K. Kobayashi, and 
H. Tsuchiura, Journal of the Physical Society of Japan 82, 
014707 (2013), http://dx.doi.org/10.7566/JPSJ.82.014707. 

[114] E. Kagawa, T. Itou, K. Miyagawa, and K. Kanoda, Phys. Rev. 
Lett. 93, 127001 (2004). 

[115] S. Sakai, G. Sangiovanni, M. Civelli, Y. Motome, K. Held, 
and M. Imada, Phys. Rev. B 85, 035102 (2012). 

[116] Y. Lu, M. M. Hoppner, O. Gunnarsson, and M. W. Haverkort, 
Phys. Rev. B 90, 085102 (2014). 

[117] L. E. Tocchio, C. Gros, X.-E. Zhang, and S. Eggert, arXiv 
preprint arXiv: 1402.3160 (2014). 

[118] H. Bakrim and C. Bourbonnais, Phys. Rev. B 90, 125119 
(2014). 

[119] PA. Igoshev, M. A. Timirgazin, A. K. Arzhnikov, and V. Y. 
Irkhin, arXiv preprint arXiv: 1409.0320 (2014). 

[120] T. Yoshikawa and M. Ogata, Phys. Rev. B 79, 144429 (2009). 


